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ABSTRACT 


The  blast  scaling  principle  by  charge  weights  may  often  produce  gross 
errors  in  the  estimation  of  overpressures  at  significant  distances  from  a 
shot  site.  The  principle  ignores  weather  conditions  and  assumes  a 
homogeneous  atmosphere.  Weather  conditions  may  cause  variations  in  the 
velocity  of  sound  which  allows  the  atmosphere  to  refract  sound  waves  and 
act  as  an  acoustic  lens.  Certain  weather  conditions  such  as  temperature 
inversions  may  allow  an  acoustic  focus  many  kilometres  away  from  a  shot 
site.  Acoustic  propagation  from  a  noise  producing  installation  has  been 
mathematically  modelled  to  assist  in  the  prediction  of  anomalous 
overpressures  by  acoustic  focusing.  -MRt-has  developed  a  computer  model, 
REFFOCUS^to  illustrate  foci  in  a  vertically  inhomogeneous,  non-turbulent, 
moving  atmosphere.  The  model  is  produced  from  sound  velocity  profiles 
derived  primarily  from  temperature  and  wind  shear  profiles.  ' 
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FOCUSING  OF  ACOUSTIC  WAVES  IN  A  NON-UNIFORM  ATMOSPHERE 


1.  INTRODUCTION 


Anomalous  overpressures,  larger  or  smaller  than  expected,  may  be  manifest  adjacent  to,  or 
at  relatively  large  distances  from,  the  boundaries  of  an  explosives  or  noise  producing 
installation.  Although  of  relatively  low  intensity,  these  focused  shock  waves  may  be 
avoided  by  ensuring  favourable  weather  conditions.  The  purpose  of  this  report  is  to 
produce  a  model  to  predict  acoustic  focused  events  as  a  function  of  the  acoustic  properties 
of  the  atmosphere. 

Normally,  assuming  a  surface  burst  (or  surface  noise  generator)  and  depending 
upon  the  amount  of  absorption  and  cratering  of  the  surface,  from  50  to  100%  of  the  energy 
generated  by  an  explosion  may  be  transmitted  radially  from  the  shot  site.  There  are 
numerous  conflicting  theoretical  and  empirical  equations  which  have  been  employed  to 
calculate  the  energy  and  overpressures  at  any  given  radius  from  the  shot  site  (Baker, 

1973).  However,  these  calculations  are  known  to  underestimate  observed  values  at  focused 
events  and  to  overestimate  observed  values  in  the  quiet  zones  in  between.  The  assumption 
of  the  continuous  divergence  of  blast  energy  at  extended  range  is  in  error  with  the 
detection  of  energy  concentrations  at  foci.  Remote  from  the  shot  site,  overpressures  are 
relatively  low  and  may  be  regarded  as  perturbations  in  the  ambient  pressure  so  that  the 
speed  of  a  front  is  approximately  acoustic  (Cox  et  al.,  1952)  and  the  front  will  obey  the  laws 
of  acoustics.  Derivations  and  calculations  in  this  report  assume  acoustic  propagation  away 
from  the  shot  site  and  are  concerned  with  the  prevailing  weather  conditions  producing  an 
acoustic  lens  within  the  layers  of  the  atmosphere. 


2.  FACTORS  AFFECTING  ACOUSTIC  FOCUSING 

2.1  Velocity  of  Sound  -  Temperature,  Wind,  Humidity  Effects 
2.1.1  Temperature  Effects 

The  equation  for  the  speed  of  sound,  v,  in  a  homogeneous  still  atmosphere  can  be  reduced  to 
(Herman  et  al.,  1983): 

1/2 


v  -  (1  +  0.26  q)  (KTy) 


(m/s) 


(1) 


where  y  -  Cp/Cv  -  ratio  of  specific  heats  at  constant  pressure  and  temperature  in  dry 

E  -  specific  gas  constant  in  dry  air  *  287.05  (J  kg-1  K-1); 

T  -  temperature  (K); 

q  -  specific  humidity  (kg/kg  of  moist  air). 


Figure  1  shows  a  variation  in  the  velocity  of  sound  due  to  atmospheric  temperature 
gradients.  An  increase  in  temperature  will  result  in  an  increased  sonic  velocity. 


UJ  Humidity  Effects 

The  contribution  to  sonic  velocity  from  humidity  is  below  1%  at  normal  temperatures  and 
may  be  neglected  (Herman  et  al). 


1L8  Wind  Shear  Effects 

Wind  velocities  may  have  a  significant  effect  on  the  velocity  of  sound  as  predicted  by  the 
above  equation.  In  a  static  atmosphere  a  sound  wave  will  be  allowed  to  travel  through  each 
medium  with  a  wave  velocity  as  prescribed  by  the  prevailing  temperature  and  humidity 
effects.  However,  in  a  non-turbulent,  moving  atmosphere  each  medium  has  a  velocity  with 
respect  to  the  ground  so  that  sound  wave  velocity  components  and  the  corresponding  wind 
velocity  components  may  be  summed  vectorially  (Fig.  1).  This  means  high  winds  may 
markedly  alter  acoustic  velocities  from  very  low  to  quite  large  sound  velocities. 

Therefore,  along  with  temperature  profile  effects,  wind  shear  effects  are  a  major 
contributor  to  the  question  of  noise  focusing  in  the  atmosphere.  The  component  of  the 
wind  velocity  which  is  of  interest  is  the  one  that  moves  along  the  path  taken  by  the 
propagating  sound  front. 


It  is  convenient  to  indicate  the  direction  of  the  energy  flow  by  a  series  of  acoustic  ray  paths 
of  equal  angular  displacement  which  radiate  out  from  the  shot  site  and  remain  at  all  times 
normal  to  the  blast  front.  In  this  fashion  an  inhomogeneous  atmosphere  may  be  indicated 
by  ray  paths  which  deviate  from  straight  radial  paths  diverging  from  the  shot  site 
(Figs  2-4).  Because  these  rays  are  traced  along  the  flow  of  energy,  then  energy  radiating 
from  the  shot  site  within  some  solid  angle  remains  within  their  bounds  -  this  region  may  be 
called  an  energy  tube  (Sach,  1961).  Although  it  is  not  anticipated  that  the  energy  within  a 
tube  will  flow  normal  to  the  bounding  rays  of  the  tube  wall,  it  is  possible  that  two  or  more 
energy  tubes  may  intersect  to  produce  a  focus.  However,  there  is  atmospheric  attenuation 
within  each  energy  tube  due  to  viscosity,  the  exchange  of  translational  and  rotational 
energy  between  colliding  molecules,  and  thermal  heating  and  conduction  effects 
particularly  near  the  shot  site  where  pressure  gradients  and  wave  velocites  are  the  greatest. 

Since  overpressure  is  related  to  the  energy  per  unit  area,  then  the  overpressure 
at  a  given  range  may  be  calculated.  If  the  energy  passes  through  a  large  area  then  the 
overpressure  is  reduced.  Conversely,  if  the  energy  passes  into  a  smaller  area  then  the 
overpressure  is  increased  to  the  extent  that  a  point  may  be  provided  with  energy  from  a 
number  of  tubes  at  varying  angles  to  produce  a  focus. 
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2 J*  Acoustic  Reflection 


The  acoustic  impedance,  Z,  of  a  material  is  a  product  of  the  velocity  of  sound,  v,  in  the 
material  and  its  density,  p. 


Z  =  pv. 


(2) 


The  fraction  of  energy  reflected  at  a  normal  to  the  boundary  between  two  media 
of  acoustic  impedances  and  Z2  is: 


(Z2  -  V 

«L  +  V2  * 


Hie  fraction  of  energy  refracted  (transmitted)  through  the  boundary  is: 


*r- 


<Z1  +  Z2> 


2  ’ 


(4) 


At  the  boundary  between  adjacent  heights  within  one  atmospheric  layer  possessing  a  linear 
velocity  profile,  or  at  the  boundary  between  two  layers  with  similar  sound  velocities  and 
acoustic  impedances  (Z^  »  Z„) ,  reflection  of  the  sound  front  normal  to  the  boundary  will 
be  a  minimum  -  e.g.  cloud  ana  meteorological  discontinuities  -  and  ER  »  0,  Erj,  *  1. 

There  is  a  maximum  of  reflection  of  the  sound  front  at  the  boundary  between 
media  of  greatly  different  acoustic  impedance  (Zi  »  Z2  or  Z1 «  Z2)  -  e.g.  air  and 
ground  -  and  Ejj  .  1,  Err,  .  0.  Note,  the  reflected  plus  refracted  energy  is  ideally  equal  to 
the  incident  energy  so  that: 


Er  ®T  =  ^ 


(5) 


depending  upon  attenuation  effects. 


2.4  Acoustic  Diffraction 

Diffraction  of  a  sound  wave  would  require  very  large  overpressures  found  in  shock  fronts 
near  the  shot  site  to  become  significant  in  the  attenuation  of  sound  wave  energy. 

However,  at  extended  ranges  the  overpressures  are  relatively  small  and  front  velocities  are 
approximately  the  speed  of  sound  (Whitham,  1957).  In  the  present  study  we  therefore 
neglect  the  effects  of  diffraction. 
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2JB  Acoustic  Refraction 

Refraction  through  the  atmosphere  contributes  to  the  focusing  of  blast  energies  at  extended 
ranges.  Section  3  details  the  equations  governing  the  path  of  sound  rays  orthogonal  to  the 
surface  of  an  acoustic  front  through  a  vertically  inhomogeneous,  stratified,  moving 
atmosphere. 

For  each  azimuth  and  within  each  atmospheric  layer  there  can  be  defined  a 
sound  velocity  profile,  K  -  dv/dy,  which  is  the  change  in  sound  velocity  with  altitude. 

Within  the  ith  atmospheric  layer  bounded  by  a  higher  altitude  sound  velocity,  vi+, ,  and  a 
lower  altitude  sound  velocity  vj,  there  will  be  a  sound  velocity  profile  Kj  where  all 
refracted  sound  rays  behave  as  follows: 


K 

Condition 

Refraction 

Figure  No. 

KjcO 

vi+l<vi 

Upward 

Figure  2 

Kj  =*  0 

vi+l  =  vi 

Not  refracted 

Figure  3 

Kj  >  0 

vi+l>vi 

Downward 

Figure  4 

Generally  a  ray  will  be  turned  toward  regions  where  the  sound  velocity  is 
lower.  According  to  a  Huygens-Kirchhoff  construction,  every  point  on  a  moving  wave  front 
of  arbitrary  shape  at  an  instant  t  is  the  source  of  secondary,  outward-propagating  spherical 
wavelets  of  variable  phase  and  of  radius  v6t .  Where  v  is  larger,  v6t  and  the  radius  of  the 
wavelet  is  larger.  The  tangents  of  the  forward-going  portions  of  these  wavelets  of 
different  radii  produce  an  envelope  which  defines  the  advancing  wave  front  at 
time  t  +  6t .  The  wave  front,  having  already  advanced  further  in  the  region  of  higher 

velocity,  will  bend  toward  the  region  of  lower  velocity. 

Ground  focusing  requires  the  presence  of  stratified  atmospheric  layers  with 
sound  velocity  profiles  such  that: 

(i)  some  sound  velocity  is  greater  at  a  higher  altitude  than  at  the  ground  -  this  is  the 
primary  condition  for  a  focused  event} 

(ii)  a  collection  of  rays  at  various  angles  are  incident  at  the  same  range,  and  arrive  at 
the  same  time  and  in-phase  (according  to  Fermat’s  principle  of  least  time)  so  that 
they  superimpose  to  produce  a  large  overpressure. 


3.  RAY  MECHANICS  OF  ATMOSPHERIC  REFRACTION 
Tie  ray  of  an  acoustic  front  has  slope,  y’: 
y'  -  -  tane 

where  0  -  angular  displacement  from  the  horizontal; 

y  -  vertical  displacement  from  the  shot  site; 


(6) 
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x  -  horizontal  displacement  from  the  shot  site. 

The  rate  of  change  in  slope,  y",  of  the  blast  front  ray  is: 


„  dy*  d(tane) 

^  *  dx  dx 


(7) 


Snell’s  law  for  two  different  media  (e.g.  two  air  layers)  describes  the 
relationship  between  velocity  and  direction  components  (Fig.  5a) 


8in(90°  -  e^)  sin(90°  -  e.+1) 


v. 

1 


vi+l 


=  c  =  constant 


where 


i  = 


i 

vi+l 


1,  2,  3 . .  n; 

is  the  number  of  atmospheric  layers; 

is  the  initial  sound  velocity  in  the  ith  layer; 

is  the  initial  sound  velocity  in  the  (i  +  l)th  layer  and  is 
equal  to  the  final  sound  velocity  in  the  ith  layer; 


Thus 


8 .  is  the  initial  angular  displacement  from  the  horizontal  in  the 
ith  layer; 

6.  -  is  the  initial  angular  displacement  from  the  horizontal  in  the 
(i  +  l)th  layer. 


cosd. 


v. 

i 


COS0  . 


i+1 


=  c 


i+1 


(8) 


Snell’s  law  may  be  iterated  through  each  layer  to  the  shot  site  so  that: 


v. 

i 


cose. 


cose. 

i 


(9) 


where:  v^  is  the  initial  sound  velocity  in  the  first  layer  (at  the  shot  site); 

91  is  the  initial  angular  displacement  from  the  horizontal  in  the  first 
layer. 


It  should  be  noted  that  Eqn.  (9)  holds  equally  well  within  individual  layers  if  the 
subscripts  are  allowed  to  identify  sound  velocities  and  angles  at  a  given  height.  Because  of 
the  relatively  large  acoustic  propagations  involved  Eqn.  (9)  assumes  that  a  blast  wave  is 
acoustic  in  behaviour  at  the  origin.  Generally,  from  figure  5b  and  Eqn  (8) 
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cose. 


(10) 


=  v.c 


and, 

sine.  =  (1  -  <v.c)2)1/2 


(11) 


If  there  exists  a  ray  which  turns  over  just  within  the  upper  boundary  of  (i  -  l)th  layer, 

thene.  =  0,  and  cose.  =  1. 
i  ’  i 

So  there  is  a  critical  angle,  ©.  .  ,  coinciding  with  the  maximum  angle  in  layer 

1  which  will  produce  a  grazing  ray  at  therop  of  the  (i  -  l)th  layer  (bottom  of  the  ith  layer), 
so  that: 


v. 

l 


cose 


lcrit 


(12a) 


Further,  any  ray  which  turns  over  in  any  layer  so  that  it  returns  to  the  ground  has  a 
velocity,  v^,  at  the  top  of  the  arc  and  must  leave  the  shot  site  at  an  angle  less  than  or  equal 
to  e  .  .  Otherwise  the  ray  will  not  turn  over  and  Cos  e.  *  1 ,  i.e.  there  is  some  sound 
velocffy  greater  than  vj  at  a  lower  altitude.  1 


So: 


v. 

l 


CO80J 


(12b) 


where  i  6lcrit  for  the  layer.  Otherwise  Gqn.  (9)  applies. 


3.1  Bay  Trace 

For  the  ith  layer,  substitute  Eqns.  (10)  and  (11)  into  Eqn.  (6), 


y’  =  tan© 


sin© 

cos© 


(1 


(vc)2) 

VC 


1/2 


(see  Figure  5b) 


=  (l/(vc)2 


1) 


1/2 


(13) 


Also  for  the  ith  layer,  substitute  Eqn.  (7),  using  partial  derivatives: 


"  _  dy*  =  dy’  d©  dy  dv 
dx  d©  dv  dx  dy 


Since  cos©  =  vc 


from  Eqn.  (10), 
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-  sinede  =  cdv 


and  hence 


de  _  c 
dv  sine  ’ 


Also  from  Eqn.  (6), 


§sL_ 

de 


sec^e 


y 


n 


=  sec  e  x 


-c  sine  dv 
sine  cv  dy 


sec  e  dv 
v  dy 


(1  +  tan^e)  dv 
v  dy 


so 


g  +  <v’)2)  dv 
v  dy 


(14) 


If  the  variation  of  the  velocity  of  sound  with  height,  dv/dy,  is  assumed  to  be 
constant  in  each  atmospheric  layer,  i.e. 


dv 

dy 


=  constant 


then 


and 


v.  ,  -  v.  =  K.  (h.  , 
i+l  i  l  i+l 


h.) 

j 


Vi+1 


=  K. 


y  +  v. 

i 


(15) 


Here  the  height,  y,  traversed  in  the  ith  atmospheric  layer  is: 


y  =  hi+1  -  hi 


(16) 
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where  hj  is  the  initial  height  at  the  origin  of  the  ray  in  the  ith  level; 
hj+1  is  the  final  height  of  the  ray  in  the  ith  level. 


By  Eqn.  (15)  the  velocity,  v,  throughout  the  layer  is: 


v  =  v.+  K.  y 


and 


(1  +  (y’)2)K. 
(K.(y  +  v./K.)) 


(1  +  (v’)2) 
y  +  v./K.  ’ 


for  K.  *  0. 


(17) 


The  case  Kj  =  0  will  be  addressed  later. 

Equation  (17)  is  a  nonlinear,  inhomogeneous  differential  equation.  A  solution 
is: 


(x  -  x.)2  +  (y  +  v./K.)2  =  r{2  (18) 

which  is  illustrated  in  figure  6a.  This  is  a  circle,  centre  =  (xj,  -  vj/Kj),  radius  =  r-,  where  xj 
and  ri  are  constants  of  integration,  x:  is  the  x  cartesian  coordinate  where  a  ray  is  at  a 
maximum  for  Kj  >  0,  or  at  a  minimum  tor  Kj  <  0  (i.e.  slope  =  0),  and  rj  is  the  radius  of 
curvature  of  the  ray. 

It  should  be  noted  that  the  centre  of  the  circle  describing  the  ray’s  path  is 
vertically  displaced  by  -v^/Kj  from  the  height  of  the  origin  of  the  ray  in  the  ith  layer,  and 
horizontally  displaced  by  x^.  If  K^  is  r  isitive,  then  x^  =  v-  tan  e .  /K.  is  +  ve,  and  -  v-  IK:  is  - 
ve;  if  K^is  negative,  then  Xj  is  -  ve  and  vj/Kj  is  +  ve.  In  cartesian  coordinates,  with  tne 
origin  being  the  point  at  which  the  ray  enters  the  atmospheric  layer  (or  the  shot  site),  a  ray 
path  with  a  positive  K  is  a  circle  with  a  centre  to  the  right  and  below  the  origin  (4th 
quadrant)  -  the  ray  is  refracted  toward  the  ground  as  shown  in  figure  6a.  A  ray  path  with  a 
negative  K  is  a  circle  with  a  centre  to  the  left  and  above  the  origin  (2nd  quadrant^  -  the  ray 
is  refracted  upward  as  shown  in  figure  6b. 

It  should  be  noted  that  there  are  theoretical  points  in  both  figures  6a  and  6b 
where  the  sound  ray  has  an  infinite  slope  and  theQvelocity  of  sound  is  zero.  This  follows 
from  Eqn.  (9)  where  the  second  angle  is  e_  =  90°  (cose„  =  0) .  Also  there  are  points 
where  the  sound  ray  turns  over  with  a  slope  of  zero  and  tne  velocity  of  sound  is  a  maximum 
-  by  Eqn.  (9)  when  ©2  =  0  (cos 9^  =  1) . 

We  now  present  examples  of  the  application  of  the  solution  to  one,  two,  three 
and  n  atmospheric  layers  where  the  K^  0. 
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(19) 


8  J..1  1  Atmospheric  Layer 
From  Eqn.  (18)  (Figure  6a  and  6b) 

<x  -  x^2  +  (y  -  hj  +  Vj^/Kj^)2  =  rj2 

where  h]_  is  the  height  of  the  origin  of  the  ray  in  the  first  layer  (the  shot  site).  Boundary 
conditions  at  the  maximum  (+  ve  K)  or  the  minimum  (-  ve  K)  are: 

x  =  x1,  y  *  hj  »  yr  so  v*  =  (yx  +  \l\^2 

■  <v2/ 

The  ray  turns  over  if  it  is  allowed  to  travel  only  in  this  layer  from  Eqns.  (12b)  and  (15): 

2 


+ 


(y  -  \  +  v1/Kl)2 


(Kj^cos  e.^) 


2  • 


(20) 


Equation  (12b)  is  true  for  any  ray  which  turns  over  in  the  ith  layer  with  an  angle 
less  than  the  critical  angle  (any  instantaneous  ray  angle  or  sound  velocity  in  the  ith  layer 
can  be  described  by  Eqn.  (8)).  Rearranging  the  circle,  Eqn.  (20)  gives 


a! _ 

ll  ~  K.  *  __ 


(KlCO801 


~2  ~  (x  "  xl) 


(21) 


where  the  square  root  is  to  be  given  the  same  sign  as  K. 


3.1.2  2  Atmospheric  Layers 

Refer  to  figure  5a  substituting  l=i,  2=i+l  with  positive  sound  velocity  profile  and  figure  7 
with  negative  sound  velocity  profiles. 

The  equation  of  the  first  circle  for  arc  AB  is  now 

2 

y 

(x  -  x.)2  +  (y  -  h.  +  v./K,)2  =  - - - 2 

(KjCos^) 

If  Kj  is  positive  then: 


xsXj 


(22) 
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where:  x^  is  now  the  horizontal  displacement  of  the  ray  in  the  first  layer,  before 
turnover; 


X^  is  the  horizontal  displacement  of  the  centre  of  the  circle  from  the  origin 
-  positive  for  >  0,  negative  for  Kj  <  0. 

Equation  (21)  becomes: 


(Kjcose^' 


-  (x  -  X1)i 


(23) 


Hie  equation  of  the  second  circle  for  arc  BC  is: 


(x  -  {xx  +  X^)2  +  (y  -  h2  +  v^)2  = 


(K2cose1) 


(24) 


If  K2  is  positive  then: 

*1  *  *  s  xx  +  X2 
and  X2  2  X2 

where:  x2  is  the  horizontal  displacement  of  the  ray,  from  the  origin  of  the  ray  in 
layer  2,  before  turnover; 

X2  is  the  horizontal  displacement  of  the  centre  of  the  circle  from  the  origin 
of  the  ray  in  layer  2  -  positive  for  K2  >  0,  negative  for  K2  <  0; 

h2  is  the  height  of  the  origin  of  the  ray  in  the  second  layer. 

So: 


y 


(x  -{Xj^ 


(25) 


8.1.3  3  Atmospheric  Layers 
Refer  to  figure  8  with  positive  sound  velocity  profiles. 

The  equation  of  the  circle  for  arc  AB  is  now: 


(x  -  Xj)2  +  (y  -  hj  +  v^/K^)2  = 


(KjCos^ 


.2  ' 
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If  Kj^  is  positive  then 


The  equation  of  the  second  circle  for  arc  BC  is: 

2 

v 

(x  -  { x1  +  Xg})2  +  (y  -  h2  +  Vg/Kg)2  = - 1 - 2  • 

If  Kg  us  positive  then:  (KgCose^ 

x^  i  x  i  x^  +  Xg 
and  Xg  i  x2 

The  equation  of  the  third  circle  for  arc  CD  is: 

2 
v 

(x  -lx±  +  Xg  +  Xg})2  +  (y  -  h3  +  Vg/Kg)2  = - - - g  (26) 

If  Kg  is  positive  then:  (KgCOse.^) 

*1  +  x2  s  x  s  xt  +  x2  +  Xg 
and  Xg  2  Xg 

Xg  is  the  horizontal  displacement  of  the  ray,  from  the  origin  of  the  ray 
in  layer  3,  before  turnover; 

Xg  is  the  horizontal  displacement  of  the  centre  of  the  circle  from  the  origin  of 
the  ray  in  layer  3  -  positive  for  Kg  >  0,  negative  for  Kg  <  0; 

hg  is  the  height  of  the  origin  of  the  ray  in  the  third  layer. 


y 

3.1.4  n  Atmospheric  Layers 

Similarly,  the  equation  of  the  nth  circle  for  the  nth  arc  is: 
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n-1 


<X  ■£Ii=lXi1  +  *n})  +  (y  -  ^  +  vn'V'  = 


(KgcosSj^)'' 


If  is  positive  then: 


n-1  n-1 

2  x{  i  t  s  [  :  x5 ]  +  X_ 
i=l  1  i=l  1  ^ 


and  2  x  and  K.»  0  (28) 

n  n  i 

where:  Xjj  is  the  horizontal  displacement  of  the  ray,  from  the  origin  of  the  ray 
in  the  nth  layer,  before  turnover; 

Xn  is  the  horizontal  displacement  of  the  centre  of  the  circle  from  the 

origin  of  the  ray  in  the  nth  layer  -  positive  for  >  0,  negative  for  Kn  <  0 

hjj  is  the  height  of  the  origin  of  the  ray  in  the  nth  layer. 


So: 


-vH 


(K  cose, )d 
n  l 


n-1 

(x  -{[  l  x.l  +  X })“ 
i=l  1 


(29) 


If  Kjj  is  positive  then 

n-1  n-1 

IXjSXSll  Xil+Xn 

i=l  1  i=l  1  ^ 

and  Xjj  2  xQ  and  Kj  -  0. 

The  square  root  term  is  to  be  given  the  same  sign  as  K. 

Equations  (19)  through  (29)  hold  for  Kj  «  0 . 

3.1.6  Layers  with  Zero  Velocity  Profiles 
If  Kj  *  0,  then  by  trigonometry,  figure  9 

V  -  For  i  .  m  <  n  (30) 

m 

Generally  for  n  layers,  if  Kj  =  0,  and  i  =  m  <  n  the  cartesian  coordinates  are: 
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(31) 


and 


h  +  x  tans  +  (h  -  h  ) 
mm  m  n  mt-i 


m-1 


(32) 


The  first  terms  in  Eqns.  (31)  and  (32)  are  calculated  using  Eqn.  (28)  up  to  the 
(m  -  l)th  layer,  and  the  second  terms  are  calculated  by  trigonometry  for  the  case  where 
Km  =  0.  (The  third  terms  may  be  calculated  using  Eqn.  (28)  for  the  (m  +  l)th  to  nth 
layers.  The  number  of  terms  is  dictated  by  the  number  of  non-adjacent  zero  velocity 
profiles. 


3.2  Range 


If  a  ray  turns  over  in  the  ith  atmospheric  layer,  the  height  traversed  by  the  ray  in  the  layer 
is  given  by  Eqn.  (16): 


y  =  h.  ,  -  h. 
J  l+l  i 


The  velocity  along  the  ray  path,  v,  after  height  y  is: 


v  =  v.  +  K.y 


(33) 


where:  v-  is  the  initial  velocity  in  the  ith  layer; 

K-  is  the  velocity  gradient  in  the  ith  layer. 

In  any  atmospheric  layer;  Eqn.  10  gives 

cos  9  -  vc. 


Now,  from  Eqn.  (6)  and  trigonometric  relations: 


and 


Then: 


tans 


sec  9  -  1 


CO80 


cose 


V(  1  +  y,S> 


(34) 
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By  geometry,  figure  5b 


y ’  =  tan©  =  ± 
and  therefore 


cv 


I - TT~ 

\(1  -  c  v  ) 


n  Xi+1 


n  yi+l 


i  M  j 

i=l  i=l  x.  i=l  y.  > 


cv 


2  2V 
(1  -  c  v  ) 


dy. 


For  the  case  when  K.  «  0 ,  substitute: 
cv  =  cos© 

cdv  =  cKjdy  =  -  sin  ©d© 

We  know  that  when 

y  =  hj,  ©  =  ©. 

and 

y  =i4+1,  ©  =  e.+r 

(Note:  The  case  for  Kj  =  0  will  be  examined  later). 

Since  the  range  of  a  ray  which  returns  to  the  same  height  is  twic 
displacements  through  all  the  atmospheric  layers,  then: 


R  =  2  l  x 
n  i=l  1 

From  Eqn.  (36) 

I 

R»-t2,lV 

-  * 2  A  «:  f 


i+1 


6i+l 


cos©sin©d© 

sin© 


cosed© . 


i-1  i  'e. 


Thus 


!0i+l 


±  2  I  -7r~  l-  sin©]  | 
i-1  **1  '0. 


±  2  l  ^7  -  sin©.+1l 


(35) 

(36) 


the  horizontal 

(37) 

(38) 

(39) 
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The  first  term  in  Eqn.  (42)  is  the  range  of  the  ray  up  to  the  mth  layer  as  calculated  by  Eqn. 
(41),  and  the  second  term  is  the  range  in  the  mth  layer  calculated  from  trigonometry 
(Fig.  9).  Equation  (41)  is  used  to  add  the  third  term  to  complete  range  calculations  until 
the  shot  site  level  is  reached.  The  number  of  terms  is  dictated  by  the  number  of  non- 
adjacent  zero  velocity  profiles. 

Equation  (42)  is  the  general  equation  for  the  range  of  the  acoustic  wave  in  an  n- 
layered  atmosphere  modified  for  K-  -  0.  It  is  used  to  complement  Eqn.  (41)  in  the 
computer  model  REFFOCUS,  described  below  in  Section  4. 
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3.3  Time  of  Arrival 


Again  the  time  of  arrival  at  the  shot  site  level  is  twice  the  time  to  reach  a  turning  point. 
According  to  Fermat's  Law  of  least  time  -  the  time  of  passage  of  a  ray  between  two  fixed 
points  is  an  extremum  (minimum)  with  respect  to  possible  paths  connecting  the  points.  The 
passage  of  sound  rays  of  variable  velocities  through  the  atmosphere  is  governed  by  Fermat’s 
law.  This  is  true  for  any  set  of  rays  (at  any  angle  ©)  which  connects  two  points  such  as  the 
shot  site  and  a  focus,  independent  of  the  path  length.  The  time  of  arrival,  tQ,  of  a  ray  back 
at  the  height  of  the  shot  site  is  given  by 


t„  =  2 


ill  I, 


li+l 


dr. 

v 


where: 


initial  height  in  the  ith  atmospheric  layer; 
final  height  in  the  (i+l)th  atmospheric  layer; 
increment  of  radial  displacement. 


The  Atmospheric  layer  thickness,  y,  is: 


(43) 


y  =  hi+i-hi. 

Hie  velocity,  v,  throughout  the  layer  is: 
v  =  vA  +  Ri  y. 

Converting  between  polar  and  rectangular  coordinates  in  2-dimensions, 


x  *  r  cos© 


y  =  r  sin© 


dr  =  - 


dx 

sin© 


dr  = 


_dy_ 

cos© 


(Note:  The  case  for  Ki  =  0  will  be  examined  later).  For  the  case  where  «  0: 


n  ^i+1 

fcn  =  2  J  J  v  .+  K  .y 
i=l  h.  • 

i 


(l/sin©)dy 
.+  K  .3 

l  iJ 


where  by  geometry  cos©  =  Vj  +  Kj  y 


dy  -  - 

sindd© 

K. 

l 

When 

y  =  h., 

©  =  ©. 
i 

we  have 

y  -  h.+1. 

6  ‘ 

(44) 
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-2  .KbS. 


1=1  i  s. 


Substitute: 


IntsecS  +  tans] 


[sees  tans  +  sec  elds 
(secs  +  tans] 

(1  +  sins]  , 
tcosS(l  +  sins)] 


Equation  (45)  thus  becomes 


n  2  r  i+1 

-strf  d“ 

1=1  l  S. 

1 


=  -  l  IntsecS  +  tanS]  I 

i=l  l  S . 

l 

n  2  sees.  +  tans. 

=  k7  ln[[secS.+1  +  tans.+1]] 


Since  this  equation  describes  rays  returning  to  the  ground,  then  the 
angle,  S.+1>  grazing  the  top  of  layer  n  is  zero,  i.e. 


Vi  =  °- 


So  the  time  of  arrival  of  a  sound  ray,  back  to  shot  site  level  is;  from  Eqn.  (48), 


-  n-1  1  [ sees .  +  tans . ] 

t  =  2 {A-  IntsecS  +  tanS  ]  +  7  ==-  In  - - - - — - — — - — r 

n  K  n  n  .  %  K.  I  sees.  ,  +  tanS.  ,] 

n  i=l  1  1+1  1+1 


Equations  (43)  through  (49)  hold  for  Kj  -  0 . 

If  Kj  =  0,  then  by  geometry,  figure  9,  for  i  =  m  <  n, 


2  v  sins  • 
m  m 


Generally  for  n  layers,  if  Kj  =  0,  and  i  =  m  <  n,  then: 
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■'m 


2^tm  +  v  sin©_ 


m 


m 


(tn" 


W}- 


(50) 


The  first  term  in  Eqn.  (50)  gives  the  time  for  the  ray  to  reach  the  mth  layer 
using  Eqn.  (49),  and  the  second  term  is  the  time  to  traverse  the  mth  layer  obtained  from 
trigonometry.  The  last  term  is  the  additional  time  for  the  ray  to  reach  shot  site  level  as 
calculated  using  Eqn.  (49).  The  number  of  terms  is  dictated  by  the  number  of  non-adjacent 
zero  velocity  profiles. 

Equation  (50)  is  the  general  equation  for  the  time  of  arrival  of  the  acoustic  wave 
in  an  n-layered  atmosphere  modified  for  =  0.  It  is  used  to  complement  Eqn.  (49)  in  the 
computer  model  REFFOCUS,  described  below  in  Section  4. 


4.  ACOUSTIC  FOCUS  PREDICTIONS 


4.1  REFFOCUS  -  A  Computerised  Aid  to  Acoustic  FOCI  Modelling 

REFFOCUS,  in  Fortran  77  code,  is  a  computerised  aid  to  predict  the  location  of  acoustic 
foci  in  the  air  by  acoustic  refraction.  Using  weather  data,  or  velocity  /height  data  it 
calculates,  tabulates  and,  if  appropriate,  plots: 

1.  The  maximum  angle,  ®lcrit»  which  will  turn  over  in  the  atmosphere. 

2.  The  range,  R,  of  rays  at  discrete  angles  from  0  =  0  to  ®lcrit  • 

3.  The  time  of  arrival,  T,  of  these  rays  from  0  =  0  to  elcrit- 

4.  The  sound  velocity/height  profiles. 

5.  The  ray  path  followed  by  each  angle  from  0=0  to  elcrit  • 

6.  Multiple  reflections  of  the  Range,  Time  of  arrival,  and  Ray  paths  for  each  angle 
from  0  =  0  to  0lcrit. 

REFFOCUS  is  based  on  a  model  which  does  not  attempt  to  use  or  model  a  real 
3-dimensional  atmosphere  but  rather  a  simplified,  horizontally  homogeneous  and  stable 
atmosphere  without  curvature  or  relief  to  produce  a  working  model  of  focused  sound  rays. 

In  accordance  with  the  working  model,  overpressures  developed  either  by  direct,  or  by 
reflected  sound  ray  convergence  have  not  been  attempted  bo  that  terrain  and  vegetation  are 
not  included  in  the  model. 


The  program  presently  assumes: 

1.  Acoustic  rather  than  shock  wave  behaviour  over  entire  range. 

2.  Surface  blasts  or  noise  generation  -  not  restricted  to  sea  level. 

3.  A  locally  flat  earth  without  curvature  or  relief. 


24 


4.  A  horizontally  homogeneous  atmosphere  where  the  weather  conditions  and  sound 
velocities  are  stable. 

5.  Rays  which  return  to  earth  at  shot  site  level  are  reflected  with  100%  of  their  energy 
from  the  flat  surface  to  return  to  earth  at  multiples  of  their  initial  range. 


REFFOCUS  will  be  extended  to  - 

1.  Calculate  approximate  overpressures  at  the  ranges  of  returned  rays. 

2.  Tabulate  all  data,  including  weather  data  and  overpressures. 

3.  Calculate  and  plot  multiple  azimuthal  foci  predictions  so  that  a  foci  contour  map  is 
generated. 

4.  Calculate  and  plot  the  foci  above  ground  which  have  observed  on  ray  traces. 

5.  Calculate  and  plot  foci  from  an  above  ground  burst. 

6.  Be  fully  transportable  to  operate  near  location  on  an  IBM  or  compatible  personal 
computer. 


REFFOCUS  is  designed  to  be  readily  used  by  an  inexperienced  person.  To 
facilitate  input  there  is  the  option  of  a  keyboard  or  data  file  input;  and  a  review  option 
allows  a  recapitulation  of  data  and  graphs. 

For  each  azimuthal  sound  velocity  profile  the  three  major  equations  used  by 
REFFOCUS  in  the  calculation  of  refracted  rays  through  n  layers  of  atmosphere  are  Eqns. 
(29),  (41),  (49). 

Equation  (29) 


=  h  - 

n 


(K  cose, )' 
n  1 


n-l 

-  (x  -{[  l  x.l  +  X })' 
i=l  1  n 


calculates  the  y  cartesian  coordinate  with  respect  to  displacement  of  the  ray  in  the  x 
direction  until  turnover  (at  x  *  xmax).  REFFOCUS  produces  an  extra  set  of  cartesian 
coordinates  by  reflection  about  xmax  so  that  the  ray’s  trajectory  is  completed  from  the  shot 
site  -  to  turnover  -  to  ground. 

Equation  (41) 


▼  H  i  ^ 

R_  -  2fc°-  tan©  +  l  [v.  tan©.  -  v.  -  1  tan©.  ,1} 
n  K  n  .  ,  K.  i  i  i  1+1 

n  i=l  i 
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calculates  the  range  for  each  ray  with  angle  ©  ■  0  to  ©  .  .  Although  x  coordinate 

calculations  are  from  the  shot  site  to  turnover,  the  ranc '  or  each  ray  is  simply  twice  x_Q_. 

lilaA 

Equation  (49) 


t 


n 


2lh~  In [sec© 
Kn  n 


+  tan©_l  + 
n 


n-1  j  sec©.  +  tan©. 

ill*i  lntsecei+1+  tan©.+11} 


calculates  the  time  of  arrival  of  each  ray  with  angle  ©  -  0  to  ©1  ...  The  total  time 

for  the  ray  to  return  to  ground  is  twice  that  of  the  time  to  tumovSr.1 

Range,  Time  of  arrival,  and  Ray  trace  coordinate  are  modified  when  Kj  =  0  as 
discussed  in  sections  3.1.5,  Eqn.  (31);  3.2,  Eqn.  (42);  3.3,  Eqn.  (50). 


4J  Typical  Velocity  Profiles 

Figures  10  to  17c  have  been  produced  by  REFFOCUS  to  illustrate  the  capabilities  of  the 
program  and  to  show  some  velocity  profiles  considered  relevant  for  a  ground  explosion. 
These  figures  are  in  addition  to  those  simple  profiles  shown  in  figures  2  to  4.  With  the 
exception  of  figures  10  and  17a,  b,  c  the  graphs  are  in  groups  of  four,  each  with  the  same 
figure  number  and  with  letters  a,  b,  c  or  d.  The  nomenclature  is  as  follows: 

a)  Sound  velocity  profile 

b)  Range  vs  initial  ray  angle 

c)  Time  of  arrival  of  sound  ray  vs  initial  ray  angle 

d)  Ray  traces 


Type  1  (Fig.  10) 


No.  of  layers  =  4 


0  m 

v^  =  340  m/s 

500  m 

v2  =  335  m/s 

750  m 

Vg  =  340  m/s 

1260  m 

v4  =  335  m/s 

Kx  =  -0.01  s"1 
K2=  0.02  b~\ 
Kg  =  -0.01  s'1 


Figure  10  shows  a  velocity  profile  where  the  sound  velocity  never  exceeds  that 
at  ground  level.  This  is  in  spite  of  a  positive  velocity  gradient  between  500  and  750 
metres.  In  this  case  REFFOCUS  will  not  generate  range,  time,  or  ray  trace  graphs  because 
no  ray  will  return  to  the  ground.  A  focus  is  not  possible. 


Type  2  (Figures  11a,  b,  c,  ed) 


No.  of  layers  -  4 


hj  -  0  m 

h2  -  500  m 
ho  m  750  m 
h4  -  1250  m 


v^  *  340  m/s 
v2  -  335  m/s 
Vg  -  345  m/s 
v4  =  340  m/s 


K,  =  -0.01  s"1 
K2=  0.04  s"; 

Kg  =  -0.01  s"1 
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Figure  11a  shows  a  velocity  profile  where  there  is  initially  a  negative  velocity 
gradient  followed  by  a  positive  gradient.  This  would  probably  arise  due  to  a  temperature 
inversion,  but  could  also  be  a  result  of  wind  shear.  Low  angle  rays  which  turn  over  before 
the  top  of  the  second  atmospheric  layer  will  be  able  to  return  to  the  ground.  This  is 
possible  because  sound  velocities  at  heights  between  625  to  750  metres  are  greater  than  the 
340  m/s  at  the  shot  site.  Higher  angle  rays  will  not  return  because  of  the  negative  gradient 
above  750  metres. 

The  maximum  or  critical  angle,  below  which  rays  may  be  allowed  to  return  to 
the  ground,  is  9.77°  in  figure  lib.  At  low  angles  there  can  be  a  large  change  in  the  final 
range  between  each  angle.  This  is  shown  as  a  large  slope.  Conversely,  near  the  critical 
angle,  e..  .  ,  the  slope  of  the  curve  is  nearly  zero  which  means  rays  of  different  angles 

are  tendlngto  converge  upon  a  focus. 

As  is  generally  the  case,  the  time  of  arrival  of  each  ray  has  a  graph  identical  to 
the  range  graph  except  for  the  units  of  measurement.  According  to  Fermat’s  law  these 
graphs  will  show  the  same  turning  points  -  maxima,  minima,  and  points  of  inflection.  By 
necessity,  these  turning  points  control  and  restrict  the  shape  of  these  two  curves  so  that 
they  are  usually  identical  -  except  when  there  is  a  single  positive  velocity  gradient  as  in 
figures  16b,  c.  Because  of  this,  it  can  be  seen  that  the  time  of  travel  of  a  ray  increases 
with  range  in  all  cases. 

Figure  lid  is  the  ray  trace.  At  ground  level  there  is  a  relative  convergence  of 
rays  between  9°  to  9.77°  at  a  range  of  about  9  km.  Rays  at  angles  greater  than  elcrit  = 
9.77°  do  not  turn  over  and  do  not  return  to  the  ground.  These  rays  are  not  shown  onany 
ray  trace  by  REFFOCUS.  Significantly,  there  may  be  a  possible  focus  between  525  to  600 
m  altitude.  This  is  not  a  ground  range,  and  as  such  is  not  considered  in  the  previous 
graphs.  The  region  where  there  are  no  returned  rays  from  0  to  9  km  is  referred  to  as  a 
quiet  zone. 


Type  3  (Figures  12a,  b,  c,  d) 
No.  of  layers  =  4 


hj  =  0  m 

v-j^  =  340  m/s 

h2  =  500  m 

v2  =  340  m/s 

=  Os' 

hg  =  750  m 

Vg  =  350  m/s 

K2  =  0.04  s~ 

h4  = 1250  m 

v4  =  350  m  / s 

Kg  =  0s_ 

There  are  two  zero  velocity  gradients  in  figure  12a  with  a  positive  gradient  in 
between.  To  allow  REFFOCUS  to  produce  range,  time  and  ray  trace  graphs  the  first 
velocity  gradient  was  made  slightly  negative  and  approaching  zero.  Any  rays  greater  than 
13.73°  in  figures  12b,  c  will  travel  into  the  third  atmospheric  layer  starting  at  750  m  and 
are  refracted  upward.  These  rays  do  not  return  to  the  ground  and  so  cannot  be  represented 
in  figures  12b,  c  and  are  not  represented  in  figure  12d.  Because  of  the  first  zero  gradient, 
rays  propagating  at  low  angles  will  radiate  radially  and  low  to  the  ground  for  nearly  infinite 
distances.  For  this  reason  REFFOCUS  allows  the  selection  of  graphical  range  and  height 
limits  to  elucidate  details  of  ray  propagation  nearer  the  shot  site.  The  rays  will  begin  to 
turn  over  in  the  second  layer  at  a  relatively  large  range  from  the  shot  site.  Figures  12b,  c 
show  that  these  very  low  angle  rays  produce  large  ranges  and  times  of  arrival  with  very 
little  chance  of  a  focus.  The  graphs  only  produce  a  zero  gradient  near  9^  .  where  a 

spread  of  rays  converge  to  produce  a  possible  focus.  This  is  more  easily  seen  in  figure  12d. 


The  larger  ranges  produced  by  the  long  radial  trajectories  of  low  angle  rays  are 
seen  in  figure  12d.  Larger  angle  rays  converge  toward  the  8  km  range.  The  ray  trace 
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graph  was  also  rescaled  to  limit  graphed  ranges  to  20  km.  Again  there  appears  to  be  a 
concentration  of  rays  in  the  atmosphere  at  a  height  of  250  to  500  m. 


Type  4  (Figures  13a,  b,  c,  d) 
No.  of  layers  =  5 


=  0  m 

Vj  =  340  m/s 

h g  =  500  m 

v2  =  342  m/s 

K,  =  0.004  s'J 

ho  =  650  m 

Vg  =  355  m/s 

K2  =  0.020  s"J 

h4  =  1150  m 

v4  =  357  m/s 

Ko  =  0.004  s-1 

hg  = 1250  m 

v5  =  357  m/s 

k4  =  0  s"1 

There  are  three  positive  velocity  gradients  in  figure  13a  up  to  1150  m,  capped  by 
a  uniform  velocity  above  it.  The  central  large  velocity  gradient  can  be  expected  to  sharply 
refract  any  rays  which  enter  the  thin  atmospheric  layer.  Figures  13b,  c  show  that  very  low 
angle  rays  are  allowed  to  propagate  small  distances  in  the  first  layer  whilst  being 
immediately  refracted  downward.  As  ray  angles  increase,  the  rays  are  allowed  to  travel 
relatively  large  distances  before  turnover  just  within  the  first  layer.  Thereafter,  ray 
angles  allow  propagation  into  the  second  layer  and,  as  the  critical  angle  for  the  second 
atmospheric  layer  is  approached  (about  15°),  the  rays  are  provided  a  more  refractive 
medium  and  are  quickly  turned  over  to  return  to  ground  at  small  ranges  and  low  time 
intervals.  These  rays  converge  toward  a  possible  focus  at  about  5.5  km  shown  by  the  flat 
region  in  figures  13b,  c  and  clearly  in  figure  13d.  Discontinuities  at  this  point  and  at  the 
boundary  of  other  layers  would  be  less  sharp  if  smaller  angle  increments  were  to  be  used  in 
the  calculation.  The  two  apparent  discontinuities  are  in  fact  the  maximum  and  minimum 
of  the  range  and  time  graphs. 

Figure  13d  shows  the  increase,  decrease,  then  increase  in  range  described  by  the 
previous  graphs.  There  is  a  rapid  increase  in  range  at  low  angles,  then  the  decrease  in 
range  as  rays  are  refracted  within  the  second  layer  and  the  final  rapid  increase  in  range  as 
the  larger  rays  near  6-  .  are  slowly  turned  over  in  the  third  layer.  Rays  of  still  greater 

angle  are  refracted  upward,  never  to  return  to  the  ground,  and  are  not  shown  in  figure  13d. 


Tjrpe  5  (Figures  14a,  b,  c,  d) 
No  of  layers  =  4 


hj  =  0  m 
h2  =  600  m 
hg  =  750  m 
h4 =  1250  m 


Vj  =  340  m/s 
v2  =  345  m/s 
Vg  =  335  m/s 
v4  =  340  m/s 


K1  =  0.01  s_1 
Kq  =  -0.04  s 
Kg  =  0.01  s_1 


There  is  only  one  atmospheric  layer  of  interest  in  figure  14a.  Any  ray  which 
has  not  turned  over  by  the  top  of  the  first  layer  will  be  refracted  upward  by  the  following 
inversion-like  layer  and  cannot  be  completely  refracted  toward  the  ground  by  the  third 
layer  of  positive  gradient.  The  third  layer  does  not  produce  sound  velocities  greater  than 
that  found  at  the  top  of  the  first  layer  and  so  does  not  contribute  toward  sound  focusing  at 
the  ground. 


Because  we  are  essentially  considering  only  a  single  atmospheric  layer  limited  to 
500  m,  figures  14b,  c  show  an  almost  linear  increase  in  range  and  time  of  arrival  with  ray 
angle  up  to  6,  . .  =  9.77°.  All  of  these  rays  are  turned  over  within  the  500  m 

ceiling.  Altn&?^n  there  is  no  focus  possible  the  sound  levels  along  the  specified  direction 
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are  higher  than  those  expected  when  acoustic  refraction  is  not  considered.  The  range  is 
limited  to  11.7  km  and  the  corresponding  time  to  34.26  s,  by  the  inversion-like  layer  above 
500  m. 


Type  6  (Figures  15a,  b,  c,  d) 
No.  of  layers  =  5 


h^  =  0  m 

v^  =  340  m/s 

h2  =  500  m 

v2  =  346  m/s 

K1=  0.012  s-1 

ho  =  650  m 

Vg  =  340  m/s 

K2  =  -0.040  s-1 

h4  = 1400  m 

v4  =  349  m/s 

Ko  =  0.012  s-1 

hg  =  1500  m 

Vg  =  349  m/s 

K4  =  0  s-1 

Hie  atmospheric  conditions  in  figure  15a  are  very  similar  to  Type  5  excepting 
the  final  velocity  gradient  is  extended  higher  so  that  there  exist  sound  velocities  greater 
than  in  lower  layers.  Rays  turning  over  are  limited  to  1400  m  where  the  last  velocity 
gradient  is  zero. 

Figs.  15b,  c  are  similar  again  to  those  of  Type  5  up  to  the  critical  angle  of  9.77° 
for  the  first  layer.  Larger  angle  rays  travel  into  the  second  layer  and  are  refracted  upward 
(Fig.  15d)  and  may  be  turned  over  in  the  third  layer  between  1250  and  1400  m.  These  rays 
travel  large  distances  in  the  second  layer  and  some  rays  with  angles  close  to  9^  .  do  turn 

over  before  some  lower  angle  rays.  There  is  a  near  focus  for  those  angles  near  eL1  .  at  a 
range  of  20  km.  cr  1 1 

The  ray  trace  below  500  m  (the  first  layer)  is  identical  to  Type  5  with  the 
addition  of  returned  rays  which  turn  over  between  1250  to  1400  m  heights.  Some  rays  have 
converged  to  a  range  of  around  20  km.  There  is  a  quiet  zone  between  10  and  20  km  since 
no  rays  have  turned  over  in  the  second  layer. 


Type  7  (Figures  16a,  b,  c,  d) 
No.  of  layers  =  2 


hjL  =  0  m  v^  =  340  m  /s 

k2  =  500  m  V2  =  345  m/s  =  0.010  s-1 


Figure  16a  has  the  single  positive  velocity  gradient  seen  in  figure  4.  Because 
the  shape  of  the  graphs  in  figures  16b,  c  are  restricted  only  by  the  one  turning  point  at  0° 
they  are  not  identical  but  there  is  still  the  asymptotic  rise  in  range  and  time  of  arrival 
near  e  .  .  Figure  16d  appears  similar  to  that  of  Type  5  because  there  is  only  one 
atmospheric  layer  of  interest  with  a  positive  velocity  gradient.  Again  these  graphs  may  be 
viewed  by  choosing  more  reasonable  range,  time  and  height  graphing  limits. 


As  in  Type  5,  sound  levels  are  enhanced  in  the  near-medium  fields  over  those 
from  direct  ray  traverse. 
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Type  8  (Figures  17a,  b,  c) 
No.  of  layers  =  5 


hl  = 

0  m 

V1 

=  340 

m/s 

1 

h2= 

500  m 

v2 

=  342 

m/s 

K1 

=  o.oo4  s'; 

h3  = 

650  m 

v3 

=  355 

m/s 

k2 

=  0.02  s'1 

h  = 

1150  m 

v4 

=  357 

m/s 

K3 

=  0.004  s'1 

bi  “ 

1250  m 

v5 

=  357 

m/s 

K4 

=  0  s-1 

These  graphs  have  been  generated  by  REFFOCUS  from  the  velocity  profile  of 
Type  4  with  one  reflection  at  the  ground  of  the  returned  rays.  Reflected  rays  may  be 
identified  in  both  the  Range  and  Time  of  Arrival  graphs  by  their  relative  amplitudes  in 
range  or  time  with  respect  to  the  unreflected  rays.  The  first  reflected  sets  of  rays  have 
twice  the  range  and  time  of  their  first  return  to  the  ground.  The  second  reflection  has 
three  times  the  range  and  time,  and  so  on.  Focii  may  experience  constructive  interference 
when  rays  of  a  different  number  of  reflections  are  incident  at  a  common  range  and  time  of 
arrival.  Figures  17a,  b  show  that  there  may  be  another  focus  derived  from  the  reflection 
of  the  first  focus.  Subsequent  foci  may  be  less  intense  compared  to  the  original  focus. 

This  may  be  shown  by  the  shorter  flat  sections  corresponding  to  reflected  rays 

of  e  m  16°  which  turn  over  just  within  the  second  layer  in  figures  17a,  b.  However,  other 

rays  may  converge  onto  the  reflected  focus,  e.g.  the  unreflected  rays  with  angles  4°  and  7°. 

Figure  17c  shows  the  ray  trace  of  both  the  unreflected  and  reflected  rays. 

Again  each  reflected  ray  may  be  identified  by  its  larger  range  than  the  corresponding 
unreflected  ray.  Also  the  arc  of  a  reflected  ray  describes  exactly  the  arc  of  an  unreflected 
ray,  beginning  at  the  return  to  ground  of  the  unreflected  ray.  A  number  of  possible  foci 
may  be  inferred  using  this  graph  and  the  range  vs  6  plot.  These  occur  at  ranges  of  5.4  and 
10.8  km,  and  at  15  km.  The  first  two  ranges  correspond  to  the  critical  angle  in  the  first 
layer.  The  ray  trace  could  be  reproduced  with  a  lower  range  limit  so  that  foci  may  be 
examined  more  closely.  Additionally  further  reflections  may  be  graphed  by  REFFOCUS. 


5.  CONCLUSION 


Installations  which  produce  noise  require  a  method  to  estimate  the  unusual  propagation  of 
noise  over  long  distances  with  anomalously  large  intensity.  The  blast  scaling  principle  by 
charge  weights  may  often  produce  gross  errors  in  the  estimation  of  overpressures  at 
significant  distances  from  a  shot  site,  the  principle  ignores  weather  conditions  and  assumes 
a  homogeneous  atmosphere.  Weather  conditions  may  cause  variations  in  the  velocity  of 
sound  which  allows  the  atmosphere  to  refract  sound  waves  and  act  as  an  acoustic  lens. 
Certain  weather  conditions  such  as  temperature  inversions  may  allow  an  acoustic  focus 
many  kilometres  away  from  a  shot  site. 

Acoustic  propagation  from  a  noise  producing  installation  has  been 
mathematically  modelled  to  assist  in  the  prediction  of  anomalous  overpressures  by  acoustic 
focusing.  Additionally,  MRL,  has  developed  a  computer  model,  REFFOCUS,  to  illustrate 
foci  in  a  vertically  inhomogeneous,  non-turbulent,  moving  atmosphere.  The  model  is 
produced  from  sound  velocity  profiles  derived  primarily  from  temperature  and  wind  shear 
profiles. 


The  results  obtained  have  shown  that  the  combined  examination  of  Range 
versus  9  and  Ray  trace  plots  for  a  given  set  of  weather  conditions  may  be  enough  to  infer 
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the  possibility  of  a  focus  by  purely  graphical  means.  This  can  be  achieved  by  considering 
the  angular  spread  of  rays  converging  down  range  from  the  shot  site.  Acoustic  foci  may  be 
predicted  by  the  collaborated  observation  of: 

1.  a  zero  slope  in  the  Range  versus  e  plot  and, 

2.  the  relative  density  of  rays  at  the  convergence  of  rays  back  to  the  ground  in  the 
Ray  trace  plot. 

Although  overpressures  are  not  estimated  in  this  model  the  contribution  of  each 
energy  tube,  containing  an  angular  spread  of  rays  incident  at  a  common  range  and  time, 
could  be  used  to  calculate  the  concentration  of  the  sound  wave  energies  at  a  focus  and  then 
related  to  overpressure.  In  this  fashion,  those  installations  which  produce  noise  or 
potentially  shocked  acoustic  foci  at  extended  ranges  may  base  the  decision  to  proceed 
according  to  REFFOCUS  predictions  using  the  prevailing  weather  data.  REFFOCUS  can 
therefore  be  used  to  minimise  noise  and  larger  overpressures  and  can  be  extended  to 
confirm  atmospheric  weather  data. 

Presently  atmospheric  data  may  be  entered  into  the  computer  model, 
REFFOCUS,  in  the  form  of  stratified  sound  velocity  data.  Alternatively,  simple 
trigonometric  relations  for  the  wind  velocity  components,  and  temperature  affected  sound 
velocity  components  has  been  included  in  the  REFFOCUS  code  to  directly  process  weather 
data.  Using  these  data  options  the  velocity  profiles,  K^,  between  each  set  of  weather 
measurements  are  calculated  and  the  full  sound  velocity  profile,  Range  versus  e,  Time  of 
arrival  of  sound  rays  versus  9,  and  Ray  traces  may  be  plotted. 

This  acoustic  focus  model  assumes  noise  generation  in  a  vertically 
inhomogeneous  atmosphere  with  stable  weather  and  sound  velocities  at  ambient  pressure. 

In  order  to  produce  a  working  acoustic  foci  model  we  have  not  considered  the  earth 
curvature,  relief,  atmospheric  composition,  blast  weights,  and  sonic  attenuation. 

This  model  of  acoustic  refraction  in  the  atmosphere,  and  REFFOCUS,  completes 
Stage  1  of  the  acoustic  focused  event  predictions.  Stage  2  will  be  used  to  experimentally 
verify  predictions  derived  from  REFFOCUS. 
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GLOSSARY  OF  SYMBOLS 


n 

i 

v 

y 

R 

T 

q 

p 

z 

Z1 

Z2 

ER 

K 

Vi 

vi+l 

t 

6t 

e 

y 

X 

/ 

y" 

9. 


Number  of  atmospheric  layers 

Subscript  i=l,2,3 . ,m  refers  to  the  ith  layer  where  msn 

Velocity  of  sound  in  air 

Cp/Cv=ratio  of  specific  heats  at  constant  pressure  and  temperature 
in  dry  air 

Specific  gas  constant  in  dry  air  =  287.05  (J  kg~*  K-1) 

Gas  temperature  (K) 

Specific  humidity  (kg/kg  of  moist  air) 

Density  of  an  acoustic  medium 
Aco  _3tic  Impedance 

Acoustic  Impedance  of  the  first  acoustic  medium 

Acoustic  Impedance  of  the  second  acoustic  medium 

The  fraction  of  energy  reflected  across  an  acoustic  boundary 

The  fraction  of  energy  transmitted  across  an  acoustic  boundary 

Velocity  profile  (=dv/dy)  of  layer  i 

Sound  velocity  at  the  foot  of  layer  i 

Sound  velocity  at  the  top  of  a  sound  ray  in  the  ith  layer  or  at 
the  foot  of  the  (i+l)th  layer 

Time 

Increment  in  time 

Angular  displacement  from  the  horizontal 
Vertical  displacement  from  the  shot  site 
Horizontal  displacement  from  the  shot  site 
Derivative  of  y,  with  respect  to  x 
Second  derivative  of  y,  with  respect  to  x 

Initial  angular  displacement  from  the  horizontal  in  the  ith  layer 


33 


B.  ...  Initial  angular  displacement  from  the  horizontal  in  the  (i+l)th  layer,  OR  final 
angular  displacement  from  the  horizontal  at  the  top  of  a  sound  ray  in  the  ith 
layer 

9.  .  Maximum  angular  displacement  from  the  horizontal  at  the  shot  site  which  will 

produce  a  grazing  ray  at  the  top  of  the  (i-l)th  layer 

hf  Initial  height  at  the  origin  of  the  ray  in  the  ith  layer 

In ,  1  Final  height  of  the  ray  in  the  ith  layer,  OR  the  initial  height  of  the  ray  in  the 

(i+l)th  layer 

yj  Vertical  displacement  of  the  ray  in  the  ith  layer 

Xj  Horizontal  displacement  of  the  ray,  from  the  origin  of  the  ray  in  the  ith  layer 

before  turnover 

Xj  Horizontal  displacement  of  the  centre  of  the  circle  from  the  origin  of  the  ray 

in  the  ith  layer 

Ru  Range  (horizontal)  of  a  ray  in  n  layer 

^  Time  of  arrival  of  a  ray  back  to  the  shot  site  height  in  n  layer 

r  Radial  displacement 

dr  Increment  of  radial  displacement 


FIGURE  3  Basic  Velocity  Profile 


FIGURE  6A  1  Atmospheric  Layer 

With  A  Positive  Sound  Velocity  Profile 


FIGURE  7  2  Atmospheric  Layers 

With  A  Negative  Sound  Velocity  Profile 


FIGURE  8  3  Atmospheric  Layers 

With  Positive  Sound  Velocity  Profiles 


FIGURE  9  An  Atmospheric  Layer  With  A  Z 


Figures  10-17  have  been  produced  using  REFFOCUS  to  illustrate  the 
capabilities  of  the  program  and  to  show  some  velocity  profiles 
considered  relevant  for  a  ground  explosion. 
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AA  a  FIGURE  11 C  TYPE  2  -  TIME  OF  ARRIVAL  vs  ELEVATION  ANGLE 


^m-)  Haight  CrrO 


Sound  Velocity  (m/s) 


Theta  1  CDeg) 


FIGURE  I 2C  TYPE  3  -  TIME  OF  ARRIVAL  vs  ELEVATION  ANGLE 
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FIGURE  158  TYPE  6  -  RANGE  vs  ELEVATION  ANGLE 
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FIGURE  168  TYPE  7  -  RANGE  vs  ELEVATION  ANGLE 
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FIGURE  17B  TYPE  8  -  TIhE  OF  ARRIVAL  vs  ELEVATION  ANGLE 
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ABSTRACT 


The  blast  scaling  principle  by  charge  weights  may  often  produce  gross  errors  in  the 
estimation  of  overpressures  at  significant  distances  from  a  shot  site.  The  principle  ignores 
weather  conditions  and  assumes  a  homogeneous  atmosphere.  Weather  conditions  may  cause 
variations  in  the  velocity  of  sound  which  allows  the  atmosphere  to  refract  sound  waves  and 
act  as  an  acoustic  lens.  Certain  weather  conditions  such  as  temperature  inversions  may 
allow  an  acoustic  focus  many  kilometres  away  from  a  shot  site.  Acoustic  propagation  from 
a  noise  producing  installation  has  been  mathematically  modelled  to  assist  in  the  prediction 
of  anomalous  overpressures  by  acoustic  focusing.  MRL  has  developed  a  computer  model, 
REFFOCUS,  to  illustrate  foci  in  a  vertically  inhomogeneous,  non- turbulent,  moving 
atmosphere.  The  model  is  produced  from  sound  velocity  profiles  derived  primarily  from 
temperature  and  wind  shear  profiles. 
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